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^ ■ ABSTRACT 

U 

^ , Planets form inside protostellar disks in a dead zone where the electrical re- 

^ \ sistivity of the gas is too high for magnetic forces to drive turbulence. We show 

that much of the dead zone nevertheless is active and flows toward the star while 
smooth, large-scale magnetic fields transfer the orbital angular momentum radi- 
ally outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of 

6 

rotational turbulence in the electrically conducting top and bottom surface layers 
of the disk, and reach the midplane by Ohmic diffusion. A toroidal component 
to the fields is produced near the midplane by the orbital shear. The process 
VO ! is similar to the magnetization of the Solar tachocline. The result is a laminar, 

^ | magnetically-driven accretion flow in the region where the planets form. 
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the dead zone gas to the magnetic fields, despite the rapid recombination of free 
charges on dust grains. Net radial magnetic fields are generated in the magneto- 



1. INTRODUCTION 

The planets of Sun-like stars form from a protostellar disk consisting of 0.01-0.1 Solar 
masses of gas together with 1% dust by mass. The removal of orbital angular momentum 
from the disk is a key to understanding the origins of the planets as it determines the rate at 
which material spirals in to accrete on the star. The magneto-rotational instability or MRI 



( IBalbus fc Hawleyill99ll ) transfers angular momentum outward and taps the free energy in 



the differential orbital rotation of the disk to drive turbulence and regenerate magnetic fields, 
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at locations where the gas is sufficiently ionized to couple strongly to the fields. However 
the minimum-mass protosolar disk, con structed by adding s ufficient hydrogen and helium to 
the planets to give Solar composition (IHayashi et al.lll985l ). is mostly too cold for thermal 
ionization. The minimum-mass disk has surface mass density S = 1700(r/AU)~ 3 / 2 g cm -2 
and temperature T = 280(r/AU)~ 1 / 2 K at radius r measured in astronomical units (AU), 
and is in vertical hydrostatic balance. Cosmic rays and stellar coronal X-rays are absorbed 
before reaching the midplane, and ionize only the surface layers. Furthermore, rapid recom- 
bination on grain surfaces leads to low abundances of the free electrons otherwise responsible 
for most of the conductivity. Consequently, much of the region where the planets form is not 
subject to the MRI. The region consists of a laminar resistive interior sandwiched between 
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is not completely quiescent, as waves propa gating from the turbule nt layers cause hydrody- 
namic stresses and low rates of accretion ([Fleming fc Stond 120031 ) . Magnetic stresses can 
occur in the interior under favorable conditions, if small grains are removed, for example by 
incorporatio n into planetesimals , and ionizing cosmic rays reach the disk unimpeded by the 
stellar wind ([Turner et al.l 120071 ). However, the resistive interior of the disk retains a weak 
coupling to the magnetic fields even if grains are present and cosmic rays are absent. The 
weak coupling is the topic of this paper. 



2. CRITERION FOR SHEAR- GENERATED MAGNETIC FIELDS 

Magnetic forces can extract angular momentum from the gas if the fields have both 
toroidal and poloidal components. We derive a criterion for the generation of toroidal fields 
by shear acting on radial fields. The Ohmic diffusion quickly eliminates any vertical magnetic 
gradients in the resistive layer, while the conducting layers above and below prevent the fields 
from escaping. The resistive layer can be approximated by an axisymmetric cylindrical flow 
in Keplerian rotation with a uniform resistivity rj. The toroidal component of the induction 
equation reduces to dB^/dt = — |QS r + r]{y 2 B ( f ) — B^/r 2 ) where B r and B^ are the radial 
and toroidal field components in cylindrical coordinates and Q is the orbital frequency. If 
also the toroidal field varies as a power-law with radius, B^ ~ r p , the field grows through 
the shear term provided r 2 Q/r] > |(p 2 — V)B^jB r . We consider the implications for the 
minimum-mass protosolar disk. The magnetic pressure is unlikely to exceed the midplane 
gas pressure, owing to the buoyancy of the fields. The magnetic pressure is a fixed fraction 
of the midplane gas pressure if the power-law index p = —13/8. The resulting criterion for 
shear-generated toroidal fields to grow to about ten times the radial seed field strength, 

v 2 J( v n) > 10, (i) 
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resembles the criterion (jSano &: Stond 120021 ) for magneto- rotational turbulence v\ z /(rjQ) > 
1. However the shear criterion equation [1] allows magnetic activity at far greater resistivities, 
since the flattened shape of the disk requires an orbital speed much greater than the sound 
speed c s , which in turn is typically greater than the vertical Alfven speed va z - The orbital 
speed is a few hundred times the vertical Alfve n speeds found at the midplane in numerical 
calculations of magneto- rotational turbulence (IMiller fc Stond 120001 ). The part of the dead 
zone where the shear criterion holds is better named the "undead zone" because it can be 
reactivated by radial fields siphoned off from the adjacent turbulent layers. In contrast with 
models where angular momentu m is removed from the top and bottom surfaces o f the disk 
by a magneto-centrifugal wind (IBlandford &: Paynd Il982l ; IWardle fc Konigll Il993l ) , here the 
torques arise from locally-generated magnetic fields. 



RESISTIVITY OF THE PROTOSOLAR DISK 



The minimum-mass protosolar disk with well-mixed micron-sized grains has an equilib- 
rium resistivity low enough for magneto-rotational turbulence only in the surface layers, at 
the radii where most of the planets formed (figured]). A deeper layer can generate smooth 
large-scale fields by shear according to equation [I] while the midplane within 5 AU of the star 
is almost completely decoupled from the fields. If the grains are removed by incorporation 
into planetary embryos, then turbulence or shear can generate magnetic fields in every part 
of the disk. 

Figured] was constructed as follows. The resistivity rj = 234v / T/x e cm 2 s _1 varies in- 
versely with the electron fraction x P = n P Jn n , where n P is the electron number density and 
n n the number density of neutrals (IBlaes fc Balbuall9941 ). We compute the electron frac- 
tion by solving a chemical network including ionization by stellar X-rays (jlgea fc Glassgold 
19991 ) and 26 Al radionuclide decay, dissociative recombination, charge exchange with metal 
atoms, and grain surface reactions. T he recombination is tre ated using the reduced reac- 
tion network with grains described by lllgner fc Nelson] (120061 ) . The fraction of the metal 
atoms free to enter the gas phase is chosen to be 1%. The resistivity depends only weakly 
on this fraction. The stellar X-ray ionizatio n rate is taken from a fit to Monte Carlo 
radiative transfer results including scattering (jlgea fc Glassgold] 119991 ). scaled to a stellar 



luminosity 2 x 10 30 erg s 1 in 5 keV thermal X-ra ys to match young Solar-mass stars 



observed in the Orion nebula ( Garmire et al] boool ). The fitted ionization rate is ( 
2.6 x 10" 15 (r/AU)- 2 (exp[-S a /8.0 g cm" 2 ] + exp[-S fe /8.0 g cm" 2 ]) s" 1 where S a and S 6 
are the mass columns lying vertically above and below the point of interest. The fit un- 
derestimates the ionization rate at mass columns < 1 g cm -2 , with no significant effect on 
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the locations of the dead and u ndead zones. T he decay of 26 Al in the dust grains yields 
an ionization rate 4 x 10~ 19 s -1 ( Stepinsk] 1992 ). A fraction 0.1% of the radioactive atoms 
is placed in the gas phase, giving a low rate of ionization even with grains absent. We ne- 
glect cosmic ray ionization owing to uncertainty about the extent to which the interstellar 
energetic particles are blocked by the wind from the young star. 



4. MHD CALCULATIONS 



We show results from two 3-D isothe rmal stratified shearing-box (IHawley et al.lll995 



Brandenburg et aDl995 : Stone et al. Il996h MHD calculations of sma ll patches of the minimum- 
mass protosolar disk using the ZEUS code (IStone fc Norman! Il992l ). The first is placed at 
5 AU and includes well-mixed 1 /im grains with dust-to-gas mass ratio 1% and a magnetic 
field with a net vertical component of 6 milligauss. The second is placed at 1 AU and has 
no grains and a net vertical field of 30 milligauss. The net magnetic fields are weak, with 
midplane ratios of gas to magnetic pressure 5 x 10 4 and 3.8 x 10 5 , respectively. The pres- 
ence of the net fields ensures that the fastest-growing mode of the linear MRI is spatially 
resolved in the calculations. To the net field in both cases is added a part with zero net 
flux that dissipates readily through reconnection. The vertical component is B sin 2irx/L, 
the toroidal component B cos2ttx/L and the radial component zero, where x is the radial 
position and L the box width. The strength B of the sinusoidal field is 0.1 and 1 Gauss, re- 
spectively. The resistivity is calculated using the ionization and chemical network described 
above except that the X-ray ionization rates are directly interpolated in the Monte Carlo 
radiative transfer results and the fit discussed above is not used. The resistivity is fixed at 
its initial equilibrium value at each height in the first run, while the resistivity is allowed 
to vary in space and time through mixing, ionization and recombin ation in the second run. 
In other respects the calculations are identical to our recent work ([Turner et al.l 120071 ). In 
particular, the vertical boundaries lie four density scale heights H = c s /Q above and below 
the midplane and are open, allowing magnetic fields to escape leaving the patch of disk with 
a net radial and toroidal magnetic flux. Among the non-ideal terms in the induction equa- 
tion, ambipolar diffusion is important outside AH. The Ohmic resistivity is greater t han the 
Hall resistivity near the midplane at 1 AU while the two are comparable at 5 AU (IWardle 
20071 ). Under these conditions the Ohmic term ha s the bigger impact on the amplitude of 
magneto- rotational turbulence (ISano fc Stondl2002l ). Our calculations therefore include only 
the Ohmic diffusion. The domain extends 2H along the radial direction and 8H along the 
toroidal direction, and is divided into 32 radial by 64 toroidal by 128 vertical zones. 



The two calculations give similar results over 130 orbits. A net radial magnetic field 
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Fig. 1. — Location of the dead zone (red) and undead zone (blue) in a cross-section view of 
the minimum-mass protosolar disk. The dead zone gas is decoupled from the magnetic fields, 
while the undead zone is sufficiently ionized for the generation of toroidal fields by shear. The 
remainder of the disk body is well-coupled to the fields and is turbulent through magnetorotational 
instability. The star lies at the origin, the disk midplane falls along the horizontal axis, and solid 
black contours show vertical mass columns of 1, 10 and 100 g cm~ 2 . The pressure in the vertical 
component of the magnetic field is set to 0.1% of the midplane gas pressure. The ionization 
is due to stellar X-rays and the decay of radioactive 26 AI. Recombination on grain surfaces is 
included in the calculation shown at top, where the dead zone extends to 5 AU in the midplane. 
The grains are 1 /im in radius and well-mixed in the gas at a 1% mass fraction. With the grains 
removed, there is no dead zone (bottom panel). 

generated in the turbulent surface layers is mixed to the undead zone where it diffuses to the 
midplane. The shear in the undead zone generates toroidal fields (figure PJ]) of the opposite 
sign, leading to magnetic stresses that transfer angular momentum outward. Lying just 
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inside the undead zone is a dissipation layer about thick that receives magnetic fields 
when turbulent motions overshoot the edge of the well-coupled gas. The resistivity is too 
high for magneto-rotational turbulence and dissipates the tangled component of the fields. 
The Ohmic diffusion time H 2 /i] at the base of the dissipation layer is about 10 orbits, so 
that the layer acts as a filter delivering to the undead zone a magnetic field that is a recent 
time-average of the large-scale field at the edge of the turbulent layer. The radial field in 
the dissipation layer and at the midplane reverses every 20-50 orbits. The midplane toroidal 
field changes at approximately the rate — % flB r due to the shea r, indicating that local shear 



generation dominates other source terms ([Turner et al.l 120071 ). The shearing reverses the 
toroidal field about 10 orbits after the radial field, so that the field lines swing around to 
trail the rotation and the magnetic stress again becomes positive (figure [3]). A related process 
occurs in the Sun, where magnetic fields generated in the con vection zone a re carried down to 
the tachocline and grow there through differential rotation (|Parkerlll993l ). The protostellar 
disk differs in that shear occurs throughout the flow, and the turbulence is shut off by the 
resistivity of the gas rather than a transition to radiative energy transport. 

The resistivity at all heights in the disk is typically unchanged as gas mixes between 
the surface layers and interior, because the recombination on the grains is fast enough to 
keep the ionization near its local equilibrium value. However wi th the grains removed , 
recombination is slow enough for ionized gas to reach the midplane (llnutsuka &: Sanoll2005l ). 
In our calculation with no grains, the undead zone shrinks after 150 years as ionized gas is 
carried down from the surface layers. The flow settles after 200 years into a new steady state 
with an overall accretion rate 20 times greater. The midplane remains magneto-rotationally 
stable (figure [3]). 



5. EMPIRICAL CONSTRAINTS 

The mass flow rate corresponding to the stresses in the MHD calculations (figure H]) is 
sufficient for the disk to accrete on the star within a few mill ion years, the lifeti me inferred 



from the disk fractions in young star clusters of different ages ( lHaisch et al.ll200ll ). The mag- 
netic contribution to the stress is several times greater than the hydrodynamic contribution 
in the dissipation layer, and at times when the fields are strongest in the midplane. The 
time-averaged midplane magnetic and hydrodynamic contributions are comparable. The 
undead zone rather than being inactive has a mean mass flow rate 4% to 61% of that in the 
turbulent layers. The midplane gas moves toward the star on average, but flows away from 
the star immediately after the reversals in the radial magnetic field when the field lines lead 
in the rotation. The root-mean-square midplane field strength found in the calculations is 



Fig. 2. — Snapshot of the toroidal magnetic field strength at 55 orbits in a resistive MHD 
calculation of a patch of the protosolar disk at 5 AU including well-mixed 1 fxm grains. The 
undead zone at center is filled with a uniform, 0.1-Gauss shear-generated toroidal magnetic field 
while patchy fields are found in the turbulent layers above and below. The star lies off-page to 
lower left and the disk midplane is horizontal through the image center. 



about 0.1 Gauss at 5 AU and 1 Gauss at 1 AU, in the range 0.1-7 Gauss inferred from the 
relict magnetization of primitive met eorites that were last m elted in the asteroid belt during 
the formation of the Solar system (ICisowski fc Hoodl Il99ll ). The agreement is consistent 
with the idea that the planet formation region is a laminar, magnetized accretion flow for 
some part of its history. 



The original dead zone picture with a fixed mass column of accreting material (iGammie 

19961 ) yields a unique mass accretion rate approximately 10~ 8 Solar masses per year, inde- 



pendent of the stellar mass. In contrast, observations of T Tauri stars show the accretion 



- 8 - 




Time (years) 



Fig. 3. — Magnetic accretion stress versus height and time in a resistive MHD calculation of 
a patch of the protosolar disk at 1 AU. The stress is proportional to the radial and toroidal 
components of the magnetic field. It is horizontally-averaged across the patch and is plotted 
using a double-logarithmic color scale, with outward angular momentum transport in red and 
yellow, and inward transport in black and blue. The calculation lasts 300 years and each panel 
shows a 100-year interval. Black curves mark the edges of the undead zone where the resistivity 
suppresses magnetorotational turbulence. Magnetic activity in the undead zone is driven by the 
diffusion of magnetic fields from the turbulent surface layers. After 150 years, the activity grows 
stronger due to an increase in magnetic coupling as ionized gas is mixed from the surface layers 
toward the midplane. 



rate increases with the stellar mass and has a spread of about two decades at Solar mass 
(IHartmann et al.ll2006l ). We suggest based on figure H] that a similar range of accretion rates 
can be produced by varying the magnetic flux and the abundance of small dust grains. 

The radial transport of material over large distances in the early S olar system is indicate d 
by the presence of crystalline silicate grains in Comet 81P/Wild 2 (IBrownlee et al.ll2006l ). 
The crystallization requires temperatures of 1000 K that were found near the Sun, while 
much of the comet is icy and formed at temperatures below 200 K in the outer Solar system. 
It remains to be seen whether grains can be carried large radial distances in the turbulent 
surface layers of the protosolar disk without sinking into the undead and dead zones. 
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Fig. 4. — Accretion rate corresponding to the total stress in the two MHD calculations located 
at 5 AU (top) and 1 AU (bottom). The accretion rate per density scale height H is plotted 
against the distance z from the midplane. The results are averaged from 20 to 120 orbits (solid 
black curves at top and bottom) and after mixing reduces the size of the undead zone, from 200 
to 300 orbits (red curve in lower panel). The vertically-integrated accretion rates are 2.2 x 10~ 8 , 
6.4 x 10 -9 and 1.5 x 10 -7 Solar masses per year respectively. Shading of the same colors marks 
the undead zone. 

6. DISCUSSION 

We have found that magnetic stresses control the dead zone evolution even with a full 
complement of micron-sized grains, in a minimum-mass protosolar disk ionized only by stellar 
X-rays and radionuclide decay. Consequently, smooth large-scale magnetic fields can occur 
in a wide range of protostellar disks and drive flows in the dead zone before and during 
planet formation as well as afterwards. 
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Models of the early evolution of the solids in protostellar disks depend crucially on the 
dynamics of the gas and will be affected by the presence of a laminar midplane accretion flow. 
The gas drag forces in turbulence cause collisions between the particle s, leading to the growth 
or de struction of solid bodies depending on the speeds involved (iDullemond fc Dominik 
20051 ). In lam inar gas, grains settle toward the midplane on timescales much less than the 
disk lifetime ( jNakagawa et al.lll98ll ). Layered accretion can cause the rapid formation of 
large bodies as particles colliding in the turbulent surface la yers grow lar ge enough to fall 
into the laminar interior where destructive collisions are rare (jCieslal 120071 ). 



Protoplanets migrate away from the locations where they first forme d, by exchanging or- 
bital angular mom entum with the disk gas through gravitational torques (jGoldreich fc Tremaine 
19791 ; I Ward! 1 199 71 ). The rapid sunward migration of the cores of Jupiter and Saturn is a dif- 



ficult y for the standard model of giant planet formation by core accretion (jHubickyj et al. 



20051 ). The migration rate de pends on the d istribution of gas near the planet, which can 
be altered by magnetic forces (ITerqueml 120031 ) . The midplane magnetic pressure in our cal- 
culations reaches a few percent of the gas pressure, sufficient for radial variations to halt 
or reverse the orbital migration. The outer edge of the undead zone is a plausible location 
for radial magnetic gradients large enough to stop the migration, and magnetic effects are a 
potential solution to the problem of the rapid loss of planetary embryos through migration. 

An important question for future studies is how the mass flow rate varies with the 
distance from the star. The fraction of the column coupled to the magnetic fields generally 
increases with the radius. Possibly the inflowing gas w ill pile up at some l ocation, leading 
to the development of a local radial pressure maximum (IKretke fc Lin! 120071 ) . If so, gas drag 
forces push solid particle s toward the maximum, pro viding favorable conditions there for the 
growth of planetesimals (jHaghighipour fc Bossll2003l ) . The further build-up of mate rial could 
lead t o self-gravitational in stability and an episode of rapid accretion on the star (iGammie 
19961 ; lArmitage et al.ll200ll ) similar to thos e associated with intervals of strong mass outflow 
from young stellar objects f|Reipurthlll989l ). 



Part of this work was carried out at the Jet Propulsion Laboratory, California Institute of 
Technology using the JPL Supercomputing Facility and with support from the JPL Research 
and Technology Development and NASA Solar Systems Origins Programs. 
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